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ABSTRACT 

Aims. Nonlinear behaviour of galactic dynamos is studied, allowing for magnetic helicity removal by the galactic fountain flow. 

Methods. A suitable advection speed is estimated, and a one-dimensional mean-field dynamo model with dynamic a-effect is explored. 
Results. It is shown that the galactic fountain flow is efficient in removing magnetic helicity from galactic discs. This alleviates the constraint 
on the galactic mean-field dynamo resulting from magnetic helicity conservation and thereby allows the mean magnetic field to saturate at a 
strength comparable to equipartition with the turbulent kinetic energy. 
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1. Introduction 

The major controversy in mean-field dynamo theory is re¬ 
lated to its nonlinear form relevant when the initial exponential 
growth of the large-scale magnetic field saturates and the field 
reaches statistical equilibrium. The core of the problem is the 
effect of the small-scale (turbulent) magnetic field on the evo¬ 
lution of the large-scale (mean) magnetic field. It has been ar¬ 
gued (Vainshtein & Cattaneo ll992l Cattaneo & Hughes fl996t 
that the Lorentz force due to the rapidly growing small-scale 
magnetic field b can make the large-scale dynamo action inef¬ 
ficient so that it produces only a negligible mean field B, with 

B ^ Rm b , where (» 1) is the magnetic Reynolds num¬ 
ber. Here and below, overbars denote averaged quantities. In 
particular, the a-effect (a key ingredient of the mean-field dy¬ 
namo) is catastrophically quenched well before the large-scale 
magnetic field can be amplified to the strength observed in as- 
trophysical objects, B ^ b. 

The suppression of the a-effect can be a consequence of the 
conservation of magnetic helicity in a medium of high electric 
conductivity (see Brandenburg & Subram anian POOS^ for a re¬ 
view). In a closed system, magnetic helicity can only evolve 
on the Ohmic time scale which is proportional to R^\ in galax¬ 
ies, this time scale by far exceeds the Hubble time. Since the 
large-scale magnetic field necessarily has non-zero helicity in 
each hemisphere through the mutual linkage of poloidal and 
toroidal fields, the dynamo also has to produce small-scale he- 
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lical magnetic fields with the opposite sign of magnetic helicity. 
Unless the small-scale magnetic field can be transported out of 
the system, it quenches the a-effect together with the mean- 
field dynamo. 

Blackman & Field J2000> first suggested that the losses 
of the small-scale magnetic helicity through the boundaries of 
the dynamo region can be essential for mean-field dynamo ac¬ 
tion. Such a helicity flux can result from the anisotropy of the 
turbulence combined with large-scale velocity shear (Vishniac 
& Cho 1200 n Subramanian & Brandenburg 120041 or the non¬ 
uniformity of the a-effect (Kleeorin et al. 120001 . The effect 
of the Vishniac-Cho flux has already been confirmed in sim¬ 
ulations, allowing the production of significant fields, B ^ b 
tBrandenburg l200-5t . 

Here we suggest another simple mechanism where the ad¬ 
vection of small-scale magnetic fields (together with the as¬ 
sociated magnetic helicity) away from the dynamo region al¬ 
lows healthy mean-field dynamo action. This effect naturally 
arises in spiral galaxies where magnetic field is generated in 
the multi-phase interstellar medium. The mean magnetic field 
is apparently produced by the motions of the warm gas (§4.3 in 
Beck et al. ll996t . which is in a state of (statistical) hydrostatic 
equilibrium with a scale height of h ^ 0.5 kpc (e.g. Korpi et 
al. ll999t . However, some of the gas is heated by supernova ex¬ 
plosions producing a hot phase whose isothermal scale height 
is 3 kpc. The hot gas leaves the galactic disc, dragging along 
the small-scale part of the interstellar magnetic field. Thus, 
the disc-halo connection in spiral galaxies represents a mecha¬ 
nism of transport of small-scale magnetic fields and small-scale 
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magnetic helicity from the dynamo active disc to the galac¬ 
tic halo. As we show here, this helps to alleviate the catas¬ 
trophic Qf-quenching under realistic parameters of the interstel¬ 
lar medium. 


2. The disc-halo connection in spiral galaxies 

The hot gas produced by supernovae cannot be confined to the 
galactic disc and therefore it is involved in systematic verti¬ 
cal motions. The gas flows to the halo where it cools and con¬ 
tracts at a height of several kiloparsecs and then falls down to 
the disc in the form of cold, relatively small clouds, forming 
what is known as the galactic fountain (Shapiro & Field 1 197^ 
Bregman 119801 . The initial vertical velocity of the hot gas is 
= 100-200 km s ' (Norman & Ikeuchi [19891 . Kahn & Brett 
[T^ Avillez & Berrv iMJH . 

The fountain flow drives gas mass flux from the disc, M - 
InR^Phfu^ through both surfaces, where / is the area covering 
factor of the hot gas, ^ 15 kpc is the radius of the galactic disc 
with vigorous supernova activity, and ph - 1-7 x 10“^^ gcm“^ 
is the hot gas density. A lower estimate of the area covering 
factor is given by the volume filling factor of the hot gas at the 
disc midplane, / = 0.2-0.3 (e.g. Korpi et al. 119991 . because 
the scale height of the galactic disc is comparable to the size of 
the hot cavities. This yields M ^ 1.5Moyr“* (cf. Norman & 
Ikeuchi [r989l who obtain M - 0.3-3 Mq yr''). 

The dynamo model discussed below refers to quantities av¬ 
eraged over scales exceeding the size of the hot cavities, and 
so it treats the multi-phase interstellar medium in an averaged 
manner. Then it is appropriate to introduce an effective foun¬ 
tain speed as the one that drives the same mass flux as above 
by advecting the diffuse interstellar gas at its mean density, 
p- 1.7 X 10“^^gcm“^ (i.e. number density of 0.1 cm“^) 

f/o =/mj—^ l-2kms^', (1) 

P 

and less for larger p (say, 0.3-0.7kms“' for the gas number 
density of 0.3 cm“^). Given the uncertainty of the magnitude of 
/ (e.g., / ^ 0.1 in the model of Norman & Ikeuchi n9891 . the 
range Uq - 0.2-2 km s“* seems to be plausible. 

The hot gas that leaves the galactic disc carries magnetic 
fields of scales smaller than the size of the hot cavities (0.1- 
1 kpc); these are mostly turbulent magnetic fields (of scales $ 
100 pc). The time scale of the removal of the small-scale mag¬ 
netic fields is of order /i/t/o - 5 x 10^ yr (with h ^ 0.5 kpc the 
scale height of the warm gas layer which hosts the mean-field 
dynamo), which is comparable to the turbulent diffusion time 
of the mean field, Irji ^ 8x10^ yr, with rji = 10^® cm^ s“' the 
turbulent magnetic diffusivity. Hence, the leakage of the small- 
scale magnetic helicity produced by the galactic fountain can 
significantly affect the mean-field dynamo. 

3. Magnetic helicity balance with fountain flow 

In the following we use, wherever appropriate, lower case char¬ 
acters to denote random, small scale quantities defined as the 
departure from the corresponding mean, e.g. a - A - A, 
j - J - J and u - U - U for the vector potential, current 


density and velocity, respectively. A gauge-invariant magnetic 
helicity density x can be defined for b as the density of corre¬ 
lated links of b using Gauss’s linkage formula (Subramanian 
& Brandenburg l2005Tl . Such a definition can be introduced for 
random magnetic fields if their correlation length is finite and 
much smaller than the system size. The helicity density a ■ b 
defined with the Coulomb gauge, V ■ a = 0, is only approxi¬ 
mately the same as x- In the presence of a helicity flux F, the 
magnetic helicity density for the small-scale field evolves as 
(Subramanian & Brandenburg l2005l 

^ + V • f = -2£ ■ B - Iripb, (2) 

at 

where £ = « x ft is the turbulent emf and rj is the Ohmic mag¬ 
netic diffusivity. (We use units where J - V y. B.) In the steady 
state, Eq. (|2li yields £> B - -jV ■ F -rjj ■ ft. Without a helicity 
flux, £> ■ B - -rjj ■ ft follows, which vanishes as rj 0 for any 
reasonable spectrum of current helicity j ■ ft so that j ■ b oc rf'‘ 
with K < 1. In the presence of helicity fluxes, however, £ ■ B 
needs not be catastrophically quenched as > 0. 

We study the effects of the simplest contribution to the flux, 

F^XU, 


which arises from the net effect of advection by upward and 
downward flows (cf. Subramanian & Brandenburg l2005l l dom¬ 
inated by the upward flow in the case of galactic fountains. The 
overall effect of the resulting helicity advection from the disc 
of the galaxy is not immediately obvious since the flow can 
also remove the mean magnetic held from the dynamo-active 
region. To reveal the net effect of the fountain flow on the dy¬ 
namo we solve the helicity equation simultaneously with the 
mean-field dynamo equation, 

^^VxiUxB + G-TjJ). ( 3 ) 

ot 

We adopt £ = aB - jjtJ (assuming isotropic turbulence) and 
take a - Ok + am (following Bouquet et al. I1976t . where ax 
represents the kinetic a-effect, am - ■ ft is the magnetic 

contribution to the a-effect, and rji - ^tu}. 

The dynamics of am is controlled by Eq. (E). We argue that 
the main contribution to am comes from the integral scale of 
the turbulence /q = Injko. Eor Kolmogorov turbulence, we 
have the following spectral scalings: ■ bt «: kb\ oc and 
Tk OC so that am oc Moreover, numerical results of 
Brandenburg & Subramanian J2005bt indicate that even j ■ ft is 
dominated by the larger scales. This justifies the estimate 

ty ~ ^ 

*^111 “ 3 * ^0 

With Bgq s pu^ and Rm - we then rewrite Eq. (|3 as (cf. 
Blackman & Brandenburg l2002l Subram anian l2002ll . 


dam 

dt 


-21],kl 


SB an 

H- 


Bi 


Rn 


- V ■ [amU) . 


(4) 


Since galactic discs are thin, it suffices to consider a one¬ 
dimensional model, retaining only the z-derivatives of the vari¬ 
ables (Ruzmaikin et al. IMH}. In terms of cylindrical coor¬ 
dinates (r, (j), z), we adopt a mean flow consisting of rotation 
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(shear) and vertical advection, with U — (0, U^, U^), where 
Uz - Uozlh within the disc. We use cck = aozih. The mean- 
field dynamo equation, in cylindrical coordinates, becomes 


dBr 

d , 


d^Br 

(5) 

dt 

dz 

[UzBz+&^) 


dB^ 

d 


, d% 


dt 

dz 

(t/,B^-£, 

l + V ‘l^oBr, 

(6) 


where q ^ -I for a flat rotation curve, oc r \ and Qq is 
the local rotation rate. We solve Eqs with the boundary 

conditions 


Br-B^-Q at z = ±h. (7) 

We adopt am = 0 at f = 0, and a symmetry condition am = 0 
at z = 0 follows automatically from the symmetry of B; the 
two conditions result in a unique solution of Eq. (|4j. (There is 
no need to be specify separately am at z = +h.) We introduce 
dimensionless numbers 


where k\ - nih. Eor the fiducial values h - 500 pc, Uq - 
0.2-2kms“', Q .0 - 25kms^'kpc^*, ao = 0.5kms“' and qt - 
10^^ cm^ s“', we have Cu — 0.1-1, Cq, ~ -2 and Ca ~ 0.8. We 
use ko/ki = fi/lo = 5 as appropriate for the galactic disc, and 
present our results in the units of 

Beq ~ 5 pG and (qik\y^ 8 x 10^ yr 

for magnetic field and time, respectively (where gas number 
density of 1 cm“^ and the turbulent velocity of 10 km s“^ has 
been used in Beq)- 

4. Results 

We compare in Eig. ^the evolution, with and without the ad- 

—2 , —2 

vective flux, of (B )/Bgq, where (B ) is the mean square of the 
large-scale magnetic field over the range |z| < h. The solutions 
are non-oscillatory and have quadrupolar parity; they grow for 
Ca > 0.2 if Cn = -2 and Cy - 0. Eor the fiducial param¬ 
eters, the e-folding time of the dynamo is about 3 x 10*^ yr for 
Cy ^ 0.3. In the absence of an advective flux, the initial growth 
of magnetic field is catastrophically quenched and the large- 
scale magnetic field decreases at about the same rate as it grew. 
The initial growth occurs while the current helicity builds up to 
cancel the kinetic a-effect. However, even a modest advective 
flux (Cy - 0.1) compensates the catastrophic quenching of the 
dynamo and the mean field energy density stays steady at about 
lO^^Bgq, which corresponds to the field strength of about 10% 
of the equipartition value. 

We illustrate in the inset of Eig.^the effect of varying the 
strength of the advection velocity on the dynamo. The steady- 
state strength of the mean field grows with Uq as the advection 
strength increases to about Cy - 0.3. Eor stronger advection, 
Cy ^ 0.5, the mean field initially grows slower but still attains 
a steady state strength slightly exceeding O.lBeq. Stronger ad¬ 
vection, Cy > 1, affects the dynamo adversely since the mean 
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Fig. 1. Evolution of the field strength at z = 0 obtained by 
solving Eqs @-(1211 with vertical advection (solid line, Cy - 
0.3) and without it (dashed line, Cy - 0), for Cn - -2, Ct, = 1 
and Bm = 10^. The dynamo is neutrally stable at Ca - 0.26 
for Cy = 0.3 and Cn = 2. The dotted curve, obtained for 
Cy <s: 1, shows that even weak advection can affect the long¬ 
term evolution of magnetic field. Eor Cy - 0, nonlinear effects 
make the a profile flatter at small |z|; this causes an oscillatory 
decay of the field. The inset shows similar results for Cy = 0.1 
(solid), 1.5 (dashed), 2 (dotted) and 3 (dash-dotted). 

field is removed too rapidly from the dynamo active region. A 
good compromise between rapid growth and large saturation 
field strength is reached for Cy ^ 0.1, which is close to the 
values expected for spiral galaxies. 

Modest advection does not noticeably affect the spatial dis¬ 
tribution of magnetic field. The profiles of B^ and B^ shown in 
Fig Efor the steady state do not differ much from solutions of 
the kinematic dynamo equations (cf. Ruzmaikin et al. 119881 . 
The corresponding profiles of a and am shown in the lower 
panel, indicate that the suppression of a in stronger near the 
disc midplane (at small z), where magnetic field is stronger. 

The steady-state strength of B, can be estimated from mag¬ 
netic helicity conservation. Averaging Eq. @, with da„ildt - 
0, over z on 0 < z < (with the mean denoted with angular 
brackets) yields 

■ B) (am) am(h)Uz(h) ^ 

Beq 277tA:„fi 

since t/j(0) = 0. Eor the adopted boundary conditions 0 
and \Bz\ \Br\, \B^\ within the disc, it is meaningful to cal¬ 
culate magnetic helicity within the disc only. Advection does 
not affect helicity conservation, which then implies {J ■ B) = 
k\{A ■ B) = -k\x- Given that £ ■ B/B^^ = (ax + am)(B/Beq)^ -i- 
am(ki/ko)^, and assuming that aK + am = ac, where ac corre¬ 
sponds to a marginally stable dynamo, Eq. 0 yields the fol¬ 
lowing estimate of magnetic field strength in the steady state: 
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Fig. 2. Plots of and Br (upper panel) and a (solid), am 
(dashed) and aK (dotted) (lower panel) versus z in the steady 
state for Cu — 0.3, Cn - -2, Ca — 1, Rm — 10^, and ki - njh. 

where D - CaCn is the dynamo number and Dc is its marginal 
value,andC = 1+Cf//(27r)+^Q/(/?m^i) = (9(1). WithA:i/A:o - 5 
and D D^., this implies B ^ 0.2Beq in rough agreement with 
the numerical solution. We note that Eq. (Unj does not apply 
for |D| » |Dc| because our arguments are only valid for weakly 
nonlinear regimes. 

5. Discussion 

The vertical advection of magnetic helicity by galactic foun¬ 
tain flow resolves straightforwardly the controversy of nonlin¬ 
ear mean-field galactic dynamos. For Cy = 0, the mean mag¬ 
netic field does initially reach a level consistent with Eq. CoJ, 
but then rapidly decays to negligible values (see Brandenburg 
& Subramanian i2005^ . The essential role of the advection is 
to provide the system with an opportunity to reach and then 
maintain the steady state CoJ within the galactic lifetime. 

These conclusions follow from Fig.^ where moderate ad¬ 
vection drastically changes the mean field levels achievable at 
t £ lO'^yr and prevents catastrophic quenching of the dy¬ 
namo. Excessive advection, however, hinders the dynamo as 
it removes the mean field from the dynamo active region. 

The steady-state strength of the mean magnetic field ob¬ 
tained in our model is of order (0.1-0.2) Beq - 0.5-1 fiG, which 
is a factor of several weaker than what is observed; we made 
no attempt here to refine the model. What is important, the 
mechanism suggested here resolves the problem of catastrophic 
quenching of the dynamo. 

The applicability of the vacuum boundary conditions Q to 
a system with gas outflow from the disc can be questionable. 
Analysis of dynamo models with outflow (Bardou et al. l200H 
suggests that reasonable changes to the boundary conditions do 
not affect the dynamo too strongly, but this aspect of our model 
should be further explored. 


We have neglected the intrinsic difference of the behaviours 
of the mean and turbulent magnetic fields near the disc sur¬ 
face. Since the horizontal size of the hot cavities is larger than 
the scale of the turbulent magnetic field but smaller than the 
scale of the mean magnetic held, the Lorentz force can resist 
the advection of the mean held more efficiently than that of the 
turbulent held. Furthermore, large-scale magnetic held loops 
drawn out by the fountain How can be detached from the par¬ 
ent magnetic lines by reconnection, so that the How will carry 
mostly small-scale Helds. Therefore, a more plausible (albeit 
less conservative) model would include advection of the small- 
scale (but not the large-scale) magnetic held. In such a model, 
the effect discussed here can be even better pronounced. 

Brandenburg et al. J1995t argue that the fountain How can 
transport the large-scale magnetic held into the halo by topo¬ 
logical pumping if the hot gas forms a percolating cluster in 
the disc and if the turbulent magnetic Reynolds number Cy in 
the fountain How exceeds ^ 20; our estimate given below is 
Cy ^ 1, and we expect that the small-scale magnetic Helds are 
removed from the disc more efficiently than the mean held. 

The idea that advection of small-scale magnetic Helds can 
help the galactic dynamo may be more robust than our partic¬ 
ular model of dynamo quenching that involves the magnetic 
a-effect. For example, if the dynamo coefficients are quenched 
due to the suppression of Lagrangian chaos by the small-scale 
magnetic Helds fKim ll99^ . their advection out of the galaxy 
will still allow the dynamo to operate efficiently. 
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